Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals.

Cycads are known to host symbiotic cyanobacteria, including Nostocales species, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition to Nostoc spp., symbiotic cyanobacteria of the genus Aulosira were identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includes Aulosira spp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.

Phenotypic heterogeneity in Streptomyces colonies.

Streptomyces are a large genus of multicellular bacteria best known for their prolific production of bioactive natural products. In addition, they play key roles in the mineralisation of insoluble resources, such as chitin and cellulose. Because of their multicellular mode of growth, colonies of interconnected hyphae extend over a large area that may experience different conditions in different parts of the colony. Here, we argue that within-colony phenotypic heterogeneity can allow colonies to simultaneously respond to divergent inputs from resources or competitors that are spatially and temporally dynamic. We discuss causal drivers of heterogeneity, including competitors, precursor availability, metabolic diversity and division of labour, that facilitate divergent phenotypes within Streptomyces colonies. We discuss the adaptive causes and consequences of within-colony heterogeneity, highlight current knowledge (gaps) and outline key questions for future studies.

Expression in CHO cells of a bacterial biosynthetic pathway producing a small non-ribosomal peptide aldehyde prevents proteolysis of recombinant proteins.

A significant problem during recombinant protein production is proteolysis. One of the most common preventive strategies is the addition of protease inhibitors, which has drawbacks, such as their short half-life and high cost, and their limited prevention of extracellular proteolysis. Actinomycetes produce the most commonly used inhibitors, which are non-ribosomal small aldehydic peptides. Previously, an unprecedented biosynthetic route involving a condensation-minus non-ribosomal peptide synthetase (NRPSs) and a tRNA utilizing enzyme (tRUE) was shown to direct the synthesis of one of these inhibitor peptides, livipeptin. Here, we show that expression of the livipeptin biosynthetic pathway encoded by the lvp genes in CHO cells resulted in the production of this metabolite with cysteine protease inhibitory activity, implying that mammalian tRNAs were recruited by the lvp system. CHO cells transiently expressing the biosynthetic pathway produced livipeptin without affecting cell growth or viability. Expression of the lvp system in CHO cells producing two model proteins, secreted alkaline phosphatase (hSeAP) and a monoclonal antibody, resulted in higher specific productivity with reduced proteolysis. We show for the first time that the expression of a bacterial biosynthetic pathway is functional in CHO cells, resulting in the efficient, low-cost synthesis of a protease inhibitor without adverse effects on CHO cells. This expands the field of metabolic engineering of mammalian cells by expressing the overwhelming diversity of actinomycetes biosynthetic pathways and opens a new option for proteolysis inhibition in bioprocess engineering.

Gut microbiomes of cycad-feeding insects tolerant to ß-methylamino-L-alanine (BMAA) are rich in siderophore biosynthesis.

Ingestion of the cycad toxins ß-methylamino-L-alanine (BMAA) and azoxyglycosides is harmful to diverse organisms. However, some insects are specialized to feed on toxin-rich cycads with apparent immunity. Some cycad-feeding insects possess a common set of gut bacteria, which might play a role in detoxifying cycad toxins. Here, we investigated the composition of gut microbiota from a worldwide sample of cycadivorous insects and characterized the biosynthetic potential of selected bacteria. Cycadivorous insects shared a core gut microbiome consisting of six bacterial taxa, mainly belonging to the Proteobacteria, which we were able to isolate. To further investigate selected taxa from diverging lineages, we performed shotgun metagenomic sequencing of co-cultured bacterial sub-communities. We characterized the biosynthetic potential of four bacteria from Serratia, Pantoea, and two different Stenotrophomonas lineages, and discovered a suite of biosynthetic gene clusters notably rich in siderophores. Siderophore semi-untargeted metabolomics revealed a broad range of chemically related yet diverse iron-chelating metabolites, including desferrioxamine B, suggesting the occurrence of an unprecedented desferrioxamine-like biosynthetic pathway that remains to be identified. These results provide a foundation for future investigations into how cycadivorous insects tolerate diets rich in azoxyglycosides, BMAA, and other cycad toxins, including a possible role for bacterial siderophores.

Streptomyces lividans 66 produces a protease inhibitor via a tRNA-utilizing enzyme interacting with a C-minus NRPS.

Small peptide aldehydes (SPAs) with protease inhibitory activity are naturally occurring compounds shown to be synthesized by non-ribosomal peptide synthetases (NRPS). SPAs are widely used in biotechnology and have been utilized as therapeutic agents. They are also physiologically relevant and have been postulated to regulate the development of their producing microorganisms. Previously, we identified an NRPS-like biosynthetic gene cluster (BGC) in Streptomyces lividans 66 that lacked a condensation (C) domain but included a tRNA-utilizing enzyme (tRUE) belonging to the leucyl/phenylalanyl (L/F) transferase family. This system was predicted to direct the synthesis of a novel SPA, which we named livipeptin. Using evolutionary genome mining approaches, here, we confirm the presence of L/F transferase tRUEs within the genomes of diverse Streptomyces and related organisms, including fusions with the anticipated C-minus NRPS-like protein. We then demonstrate genetic functional cooperation between the identified L/F-transferase divergent tRUE homolog with the C-minus NRPS, leading to the synthesis of a metabolic fraction with protease inhibitory activity. Semisynthetic assays in the presence of RNAse revealed that the productive interaction between the tRUE and the C-minus NRPS enzymes is indeed tRNA dependent. We expect our findings to boost the discovery of SPAs, as well as the development of protease-mediated biotechnologies, by exploiting the uncovered genetic basis for synthesizing putative acetyl-leu/phe-arginine protease inhibitors. Furthermore, these results will facilitate the purification and structural elucidation of livipeptin, which has proven difficult to chemically characterize. SIGNIFICANCE: The discovery of natural products biosynthetic genes marks a significant advancement in our understanding of these metabolites, for example of their evolution, activity, and biosynthesis, but also opens biotechnological opportunities and knowledge to advance genome mining approaches. We made this possible by uncovering a new biosynthetic pathway in Streptomyces lividans 66 shown to direct the synthesis of a strong protease inhibitor, termed livipeptin, following unprecedented biosynthetic rules and genes. Thus, by shedding light on the genetic mechanisms predicted to govern the production of acetyl-leu/phe-arginine protease inhibitors, including the elusive livipeptin, this study enables novel protease-mediated biotechnologies as well as approaches for discovering protease inhibitors from genome data.

On the evolution of natural product biosynthesis.

Natural products are the raw material for drug discovery programmes. Bioactive natural products are used extensively in medicine and agriculture and have found utility as antibiotics, immunosuppressives, anti-cancer drugs and anthelminthics. Remarkably, the natural role and what mechanisms drive evolution of these molecules is relatively poorly understood. The exponential increase in genome and chemical data in recent years, coupled with technical advances in bioinformatics and genetics have enabled progress to be made in understanding the evolution of biosynthetic gene clusters and the products of their enzymatic machinery. Here we discuss the diversity of natural products, incorporating the mechanisms that govern evolution of metabolic pathways and how this can be applied to biosynthetic gene clusters. We build on the nomenclature of natural products in terms of primary, integrated, secondary and specialised metabolism and place this within an ecology-evolutionary-developmental biology framework. This eco-evo-devo framework we believe will help to clarify the nature and use of the term specialised metabolites in the future.

Biosynthetic novelty index reveals the metabolic potential of rare actinobacteria isolated from highly oligotrophic sediments.

Calculations predict that testing of 5 000-10 000 molecules and >1 billion US dollars (£0.8 billion, £1=$1.2) are required for one single drug to come to the market. A solution to this problem is to establish more efficient protocols that reduce the high rate of re-isolation and continuous rediscovery of natural products during early stages of the drug development process. The study of 'rare actinobacteria' has emerged as a possible approach for increasing the discovery rate of drug leads from natural sources. Here, we define a simple genomic metric, defined as biosynthetic novelty index (BiNI), that can be used to rapidly rank strains according to the novelty of the subset of encoding biosynthetic clusters. By comparing a subset of high-quality genomes from strains of different taxonomic and ecological backgrounds, we used the BiNI score to support the notion that rare actinobacteria encode more biosynthetic gene cluster (BGC) novelty. In addition, we present the isolation and genomic characterization, focused on specialized metabolites and phenotypic screening, of two isolates belonging to genera Lentzea and Actinokineospora from a highly oligotrophic environment. Our results show that both strains harbour a unique subset of BGCs compared to other members of the genera Lentzea and Actinokineospora. These BGCs are responsible for potent antimicrobial and cytotoxic bioactivity. The experimental data and analysis presented in this study contribute to the knowledge of genome mining analysis in rare actinobacteria and, most importantly, can serve to direct sampling efforts to accelerate early stages of the drug discovery pipeline.

MIBiG 3.0: a community-driven effort to annotate experimentally validated biosynthetic gene clusters.

With an ever-increasing amount of (meta)genomic data being deposited in sequence databases, (meta)genome mining for natural product biosynthetic pathways occupies a critical role in the discovery of novel pharmaceutical drugs, crop protection agents and biomaterials. The genes that encode these pathways are often organised into biosynthetic gene clusters (BGCs). In 2015, we defined the Minimum Information about a Biosynthetic Gene cluster (MIBiG): a standardised data format that describes the minimally required information to uniquely characterise a BGC. We simultaneously constructed an accompanying online database of BGCs, which has since been widely used by the community as a reference dataset for BGCs and was expanded to 2021 entries in 2019 (MIBiG 2.0). Here, we describe MIBiG 3.0, a database update comprising large-scale validation and re-annotation of existing entries and 661 new entries. Particular attention was paid to the annotation of compound structures and biological activities, as well as protein domain selectivities. Together, these new features keep the database up-to-date, and will provide new opportunities for the scientific community to use its freely available data, e.g. for the training of new machine learning models to predict sequence-structure-function relationships for diverse natural products. MIBiG 3.0 is accessible online at https://mibig.secondarymetabolites.org/.

Evolutionary Genome Mining for the Discovery and Engineering of Natural Product Biosynthesis.

Genome mining has become an invaluable tool in natural products research to quickly identify and characterize the biosynthetic pathways that assemble secondary or specialized metabolites. Recently, evolutionary principles have been incorporated into genome mining strategies in an effort to better assess and prioritize novelty and understand their chemical diversification for engineering purposes. Here, we provide an introduction to the principles underlying evolutionary genome mining, including bioinformatic strategies and natural product biosynthetic databases. We introduce workflows for traditional genome mining, focusing on the popular pipeline antiSMASH, and methods to predict enzyme substrate specificity from genomic information. We then provide an in-depth discussion of evolutionary genome mining workflows, including EvoMining, CORASON, ARTS, and others, as adopted by our group for the discovery and prioritization of natural products biosynthetic gene clusters and their products.

Microbial Diversity in Cultivated and Feral Vanilla Vanilla planifolia Orchids Affected by Stem and Rot Disease.

The worldwide production of vanilla, a native orchid from Mexico, is greatly affected by stem and root rot disease (SRD), typically associated with Fusarium oxysporum fungi. We hypothesized that the presence of Fusarium species in vanilla is not sufficient for the plant to express symptoms of the disease. We described the taxonomic composition of endophytic microbiomes in symptomatic and asymptomatic vanilla plants using 16S and ITS rDNA metabarcoding, and ITS Sanger sequences generated from fungal isolates. We compared the bacterial and fungal diversity in vanilla plants from a long-term plantation, and from feral plants found near abandoned plantations that did not present SRD symptoms. No significant differences were found in the species richness of the bacterial and fungal microbiome among feral, or asymptomatic and symptomatic cultivated vanilla. However, significant differences were detected in both fungal and bacterial diversity from different organs in the same plant, with roots being more diverse than stems. We found that Proteobacteria and Actinobacteria, as well as the fungal families Nectriaceae and Xylariaceae, constitute the core of the vanilla microbiome that inhabits the root and stem of both cultivated and feral plants. Our work provides information on the microbial diversity associated to root and stem rot in vanilla and lays the groundwork for a better understanding of the role of the microbiome in vanilla fungal diseases.

Phylogenomics and population genomics of SARS-CoV-2 in Mexico during the pre-vaccination stage reveals variants of interest B.1.1.28.4 and B.1.1.222 or B.1.1.519 and the nucleocapsid mutation S194L associated with symptoms.

Understanding the evolution of the SARS-CoV-2 virus in various regions of the world during the Covid-19 pandemic is essential to help mitigate the effects of this devastating disease. We describe the phylogenomic and population genetic patterns of the virus in Mexico during the pre-vaccination stage, including asymptomatic carriers. A real-time quantitative PCR screening and phylogenomic reconstructions directed at sequence/structure analysis of the spike glycoprotein revealed mutation of concern E484K in genomes from central Mexico, in addition to the nationwide prevalence of the imported variant 20C/S:452R (B.1.427/9). Overall, the detected variants in Mexico show spike protein mutations in the N-terminal domain (i.e. R190M), in the receptor-binding motif (i.e. T478K, E484K), within the S1-S2 subdomains (i.e. P681R/H, T732A), and at the basis of the protein, V1176F, raising concerns about the lack of phenotypic and clinical data available for the variants of interest we postulate: 20B/478K.V1 (B.1.1.222 or B.1.1.519) and 20B/P.4 (B.1.1.28.4). Moreover, the population patterns of single nucleotide variants from symptomatic and asymptomatic carriers obtained with a self-sampling scheme confirmed the presence of several fixed variants, and differences in allelic frequencies among localities. We identified the mutation N:S194L of the nucleocapsid protein associated with symptomatic patients. Phylogenetically, this mutation is frequent in Mexican sub-clades. Our results highlight the dual and complementary role of spike and nucleocapsid proteins in adaptive evolution of SARS-CoV-2 to their hosts and provide a baseline for specific follow-up of mutations of concern during the vaccination stage.

The confluence of big data and evolutionary genome mining for the discovery of natural products.

This review covers literature between 2003-2021The development and application of genome mining tools has given rise to ever-growing genetic and chemical databases and propelled natural products research into the modern age of Big Data. Likewise, an explosion of evolutionary studies has unveiled genetic patterns of natural products biosynthesis and function that support Darwin's theory of natural selection and other theories of adaptation and diversification. In this review, we aim to highlight how Big Data and evolutionary thinking converge in the study of natural products, and how this has led to an emerging sub-discipline of evolutionary genome mining of natural products. First, we outline general principles to best utilize Big Data in natural products research, addressing key considerations needed to provide evolutionary context. We then highlight successful examples where Big Data and evolutionary analyses have been combined to provide bioinformatic resources and tools for the discovery of novel natural products and their biosynthetic enzymes. Rather than an exhaustive list of evolution-driven discoveries, we highlight examples where Big Data and evolutionary thinking have been embraced for the evolutionary genome mining of natural products. After reviewing the nascent history of this sub-discipline, we discuss the challenges and opportunities of genomic and metabolomic tools with evolutionary foundations and/or implications and provide a future outlook for this emerging and exciting field of natural product research.

Catalytically impaired TrpA subunit of tryptophan synthase from Chlamydia trachomatis is an allosteric regulator of TrpB.

Intracellular growth and pathogenesis of Chlamydia species is controlled by the availability of tryptophan, yet the complete biosynthetic pathway for l-Trp is absent among members of the genus. Some representatives, however, preserve genes encoding tryptophan synthase, TrpAB - a bifunctional enzyme catalyzing the last two steps in l-Trp synthesis. TrpA (subunit α) converts indole-3-glycerol phosphate into indole and glyceraldehyde-3-phosphate (α reaction). The former compound is subsequently used by TrpB (subunit ß) to produce l-Trp in the presence of l-Ser and a pyridoxal 5'-phosphate cofactor (ß reaction). Previous studies have indicated that in Chlamydia, TrpA has lost its catalytic activity yet remains associated with TrpB to support the ß reaction. Here, we provide detailed analysis of the TrpAB from C. trachomatis D/UW-3/CX, confirming that accumulation of mutations in the active site of TrpA renders it enzymatically inactive, despite the conservation of the catalytic residues. We also show that TrpA remains a functional component of the TrpAB complex, increasing the activity of TrpB by four-fold. The side chain of non-conserved ßArg267 functions as cation effector, potentially rendering the enzyme less susceptible to the solvent ion composition. The observed structural and functional changes detected herein were placed in a broader evolutionary and genomic context, allowing identification of these mutations in relation to their trp gene contexts in which they occur. Moreover, in agreement with the in vitro data, partial relaxation of purifying selection for TrpA, but not for TrpB, was detected, reinforcing a partial loss of TrpA functions during the course of evolution.

ActDES - a curated Actinobacterial Database for Evolutionary Studies.

Actinobacteria is a large and diverse phylum of bacteria that contains medically and ecologically relevant organisms. Many members are valuable sources of bioactive natural products and chemical precursors that are exploited in the clinic and made using the enzyme pathways encoded in their complex genomes. Whilst the number of sequenced genomes has increased rapidly in the last 20 years, the large size, complexity and high G+C content of many actinobacterial genomes means that the sequences remain incomplete and consist of large numbers of contigs with poor annotation, which hinders large-scale comparative genomic and evolutionary studies. To enable greater understanding and exploitation of actinobacterial genomes, specialized genomic databases must be linked to high-quality genome sequences. Here, we provide a curated database of 612 high-quality actinobacterial genomes from 80 genera, chosen to represent a broad phylogenetic group with equivalent genome re-annotation. Utilizing this database will provide researchers with a framework for evolutionary and metabolic studies, to enable a foundation for genome and metabolic engineering, to facilitate discovery of novel bioactive therapeutics and studies on gene family evolution. This article contains data hosted by Microreact.

Desert Environments Facilitate Unique Evolution of Biosynthetic Potential in Streptomyces.

Searching for new bioactive metabolites from the bacterial genus Streptomyces is a challenging task. Combined genomic tools and metabolomic screening of Streptomyces spp. native to extreme environments could be a promising strategy to discover novel compounds. While Streptomyces of desertic origin have been proposed as a source of new metabolites, their genome mining, phylogenetic analysis, and metabolite profiles to date are scarcely documented. Here, we hypothesized that Streptomyces species of desert environments have evolved with unique biosynthetic potential. To test this, along with an extensive characterization of biosynthetic potential of a desert isolate Streptomyces sp. SAJ15, we profiled phylogenetic relationships among the closest and previously reported Streptomyces of desert origin. Results revealed that Streptomyces strains of desert origin are closer to each other and relatively distinct from Streptomyces of other environments. The draft genome of strain SAJ15 was 8.2 Mb in size, which had 6972 predicted genes including 3097 genes encoding hypothetical proteins. Successive genome mining and phylogenetic analysis revealed the presence of putative novel biosynthetic gene clusters (BGCs) with low incidence in another Streptomyces. In addition, high-resolution metabolite profiling indicated the production of arylpolyene, terpenoid, and macrolide compounds in an optimized medium by strain SAJ15. The relative abundance of different BGCs in arid Streptomyces differed from the non-arid counterparts. Collectively, the results suggested a distinct evolution of desert Streptomyces with a unique biosynthetic potential.

Evolutionary dynamics of natural product biosynthesis in bacteria.

Covering: 2008 up to 2019The forces of biochemical adaptive evolution operate at the level of genes, manifesting in complex phenotypes and the global biodiversity of proteins and metabolites. While evolutionary histories have been deciphered for some other complex traits, the origins of natural product biosynthesis largely remain a mystery. This fundamental knowledge gap is surprising given the many decades of research probing the genetic, chemical, and biophysical mechanisms of bacterial natural product biosynthesis. Recently, evolutionary thinking has begun to permeate this otherwise mechanistically dominated field. Natural products are now sometimes referred to as 'specialized' rather than 'secondary' metabolites, reinforcing the importance of their biological and ecological functions. Here, we review known evolutionary mechanisms underlying the overwhelming chemical diversity of bacterial secondary metabolism, focusing on enzyme promiscuity and the evolution of enzymatic domains that enable metabolic traits. We discuss the mechanisms that drive the assembly of natural product biosynthetic gene clusters and propose formal definitions for 'specialized' and 'secondary' metabolism. We further explore how biosynthetic gene clusters evolve to synthesize related molecular species, and in turn how the biological and ecological roles that emerge from metabolic diversity are acted on by selection. Finally, we reconcile chemical, functional, and genetic data into an evolutionary model, the dynamic chemical matrix evolutionary hypothesis, in which the relationships between chemical distance, biomolecular activity, and relative fitness shape adaptive landscapes.

A computational framework to explore large-scale biosynthetic diversity.

Genome mining has become a key technology to exploit natural product diversity. Although initially performed on a single-genome basis, the process is now being scaled up to mine entire genera, strain collections and microbiomes. However, no bioinformatic framework is currently available for effectively analyzing datasets of this size and complexity. In the present study, a streamlined computational workflow is provided, consisting of two new software tools: the 'biosynthetic gene similarity clustering and prospecting engine' (BiG-SCAPE), which facilitates fast and interactive sequence similarity network analysis of biosynthetic gene clusters and gene cluster families; and the 'core analysis of syntenic orthologues to prioritize natural product gene clusters' (CORASON), which elucidates phylogenetic relationships within and across these families. BiG-SCAPE is validated by correlating its output to metabolomic data across 363 actinobacterial strains and the discovery potential of CORASON is demonstrated by comprehensively mapping biosynthetic diversity across a range of detoxin/rimosamide-related gene cluster families, culminating in the characterization of seven detoxin analogues.

Comparative Genomics and Metabolomics Analyses of Clavulanic Acid-Producing Streptomyces Species Provides Insight Into Specialized Metabolism.

Clavulanic acid is a bacterial specialized metabolite, which inhibits certain serine ß-lactamases, enzymes that inactivate ß-lactam antibiotics to confer resistance. Due to this activity, clavulanic acid is widely used in combination with penicillin and cephalosporin (ß-lactam) antibiotics to treat infections caused by ß-lactamase-producing bacteria. Clavulanic acid is industrially produced by fermenting Streptomyces clavuligerus, as large-scale chemical synthesis is not commercially feasible. Other than S. clavuligerus, Streptomyces jumonjinensis and Streptomyces katsurahamanus also produce clavulanic acid along with cephamycin C, but information regarding their genome sequences is not available. In addition, the Streptomyces contain many biosynthetic gene clusters thought to be "cryptic," as the specialized metabolites produced by them are not known. Therefore, we sequenced the genomes of S. jumonjinensis and S. katsurahamanus, and examined their metabolomes using untargeted mass spectrometry along with S. clavuligerus for comparison. We analyzed the biosynthetic gene cluster content of the three species to correlate their biosynthetic capacities, by matching them with the specialized metabolites detected in the current study. It was recently reported that S. clavuligerus can produce the plant-associated metabolite naringenin, and we describe more examples of such specialized metabolites in extracts from the three Streptomyces species. Detailed comparisons of the biosynthetic gene clusters involved in clavulanic acid (and cephamycin C) production were also performed, and based on our analyses, we propose the core set of genes responsible for producing this medicinally important metabolite.

Unlocking a high bacterial diversity in the coralloid root microbiome from the cycad genus Dioon.

Cycads are among the few plants that have developed specialized roots to host nitrogen-fixing bacteria. We describe the bacterial diversity of the coralloid roots from seven Dioon species and their surrounding rhizosphere and soil. Using 16S rRNA gene amplicon sequencing, we found that all coralloid roots are inhabited by a broad diversity of bacterial groups, including cyanobacteria and Rhizobiales among the most abundant groups. The diversity and composition of the endophytes are similar in the six Mexican species of Dioon that we evaluated, suggesting a recent divergence of Dioon populations and/or similar plant-driven restrictions in maintaining the coralloid root microbiome. Botanical garden samples and natural populations have a similar taxonomic composition, although the beta diversity differed between these populations. The rhizosphere surrounding the coralloid root serves as a reservoir and source of mostly diazotroph and plant growth-promoting groups that colonize the coralloid endosphere. In the case of cyanobacteria, the endosphere is enriched with Nostoc spp and Calothrix spp that are closely related to previously reported symbiont genera in cycads and other early divergent plants. The data reported here provide an in-depth taxonomic characterization of the bacterial community associated with coralloid root microbiome. The functional aspects of the endophytes, their biological interactions, and their evolutionary history are the next research step in this recently discovered diversity within the cycad coralloid root microbiome.